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Abstract 

Andreev reflection is a smart tool to investigate the spin polarisation P of the current through 
point contacts between a superconductor and a ferromagnet. We compare different models to 
extract P from experimental data and investigate the dependence of P on different contact param- 
eters. 
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Following the pioneering work of Igor Yanson and his group at the Institute for Low Tem- 
perature Physics and Engineering of the Ukrainian Academy of Sciences (ILTPE NASU) on 
tiny metallic contacts between two metal electrodes, 1 point-contact spectroscopy (PCS) has 
become a powerful method to study the interactions of ballistic electrons with other exci- 
tations in metals. 2 The interpretation of the observed characteristics in point-contact (PC) 
spectra is usually difficult. These difficulties are frequently inherent in the fabrication of 
point contacts. In many cases contacts are made by the needle-anvil or shear technique in 
which two sharpened metal pieces are brought into a gentle touch until a conductive contact 
is formed. Those contacts are microscopically not well-defined with respect to contact size, 
geometry, and structure of the metallic nanobridge, and with respect to the local electronic 
parameters such as the mean free path in the immediate contact region. The only control 
parameter is the contact resistance, and hence, it is challenging to identify the relevant 
transport regime free of doubt. Usually Sharvin's 3 or Wexler's 4 formulae for the ballistic 
and diffusive transport regime, respectively, are used to infer a PC size estimate from the 
measured PC resistance. Only recently, 5 direct scanning electron microscopy (SEM) mea- 
surements of the nanocontact size of nanostructured point contacts allowed for the first 
time a direct comparison with theoretical models for contact-size estimates of heterocon- 
tacts. The semiclassical models yield reasonable values for the PC radius a as long as the 
correct transport regime is determined by taking into account the local transport parameters 
of the individual contact. Of course, this requires a careful characterisation of the samples 
with respect to the local resistivity and the local mean free path. 

Among the rich variety of solid-state problems investigated by point-contact spectroscopy 
the study of superconductor-metal contacts contributes a significant portion. Nowadays 
point-contact spectroscopy is an important tool to explore the symmetry and nodal struc- 
ture of the energy gap A of conventional and unconventional superconductors. 6 When the 
temperature is lowered below the superconducting transition temperature T c of the super- 
conducting electrode of a superconductor (S)/ normal metal (N) point contact Andreev 
reflection 7 of charge carriers at the S/N interface occurs. Andreev reflection leads to min- 
ima at V ~ ±A/e in the differential resistance dV/dl as a function of applied bias V, i.e. 
maxima in the corresponding conductance curves G(V) = (dl /dV)(V), and thus allows 
determination of the gap size, also while varying temperature and magnetic field, respec- 
tively. 8 ' 9 
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A new pitch came into the field when Andreev reflection was used to extract the spin 
polarisation P of the current through superconductor/ferromagnet (F) point contacts. 10,11 
Knowledge of the spin polarisation of possible materials for spin-electronic devices is a key 
issue for spintronics. 12 An efficient spin injection is of central importance for utilizing the spin 
degree of freedom as a new functionality in spin-electronic devices. The spin polarisation 
P of ferro magnets can be measured by various techniques including photoemission, 13 spin- 
dependent tunnelling, 14 and point-contact Andreev reflection (PCAR). 15 For a quantitative 
analysis of the results, however, one has to be aware of the different nature of the quantities 
measured by each technique. The spin polarisation, defined by the difference of spin-up and 
spin-down density-of-states, is typically measured with spin-polarised photoemission while 
the spin polarisation of the transport current is obtained, e. g. in PCAR experiments, 16 which 
in turn is distinctly different from the spin polarisation of the density-of-states resulting from 
tunnelling experiments. 14 An issue of considerable importance is how the spin polarisation 
obtained by Andreev reflection is related to the ferromagnet's bulk spin polarisation. 17 

Already a variety of materials have been investigated including the ferromagnetic elements 
Fe, Co, and Ni and several alloys mainly with Al, Nb, or Pb as a superconducting counter 
electrode. 10 > n > 18 - 23 However, different models 10,11,24-26 describing the transport through S/F 
interfaces yielded varying values for P, also depending on the contact fabrication and the 
transport regime, 27 an issue that is not yet understood in detail. 28 In the following, we want 
to review the main ideas of two most prominent models shortly and compare the results of 
both to the same set of experimental data obtained on nanostructured Al/Fe contacts. 23 

The theoretical analysis of most S/F point-contact experiments has been carried out in 
the spirit of the Blonder-Tinkham-Klapwijk (BTK) theory 9 for Andreev reflection at an 
interface between N and classical S with spin-singlet pairing. This is the coherent process 
by which an electron from N enters S and a hole of opposite spin is retro-reflected, creating 
a spin-singlet Cooper pair in S. Possible ordinary reflection at the S/F interface barrier is 
parametrised by a phenomenological parameter, the barrier strength Z. The sensitivity of 
the Andreev process to the spin of the carriers originates from the conservation of the spin 
direction at the interface. Consequently, when there is an imbalance in the number of spin- 
up and spin-down electrons at the Fermi level, as it is the case in the spin-polarised situation 
of a ferromagnetic metal, this leads to a reduction of the Andreev reflection probability. 15 
Andreev reflection is limited by the minority carriers of the metal. 
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In the simpliest approach applied for the analysis of several experiments, 10,18-21,24 the total 
current / through the constriction is decomposed into a fully unpolarised part (1 — P')I U for 
which Andreev reflection is allowed and into a fully polarised part P'l p for which Andreev 
reflection is zero, 

I=(1-P')I U + P'I P . 

The weighting factor P 1 determines the spin polarisation of the ferromagnet. In the following 
we will refer to this model as the dispartment model. Consequentially, the conductance Gsf 
is also decomposed into two parts: 

G SF = (1 - P')G U {V) + P'G P (V) 

with 

G U , P (V) = J™ df{E ~ v V ' T) [l + A U , P (E, Z) - B U , P (E, Z)]dE 

where G u denotes the conductance, A u the Andreev reflection probability, and B u the nor- 
mal reflection probability of the fully unpolarised channel, and G p , A p , and B p denote the 
corresponding quantities of the fully polarised channel. Both contributions are derived in 
the BTK formalism and following expressions for the zero-temperature conductances G u and 
G n are obtained: 25 





\E\<A \E\>A 


G u \t=o(E) 
G P \ T=0 (E) 


2(l+/3 2 ) 2/3 


/3 2 + (l+2Z 2 ) 2 1+/3+2Z 2 
U (l+/3) 2 +4Z 2 



with p = E/^/\A 2 -E 2 \. 

Despite the attractive simplicity of the BTK formalism it has been shown 17 ' 27 that ap- 
plication of the BTK formalism (even in its generalized form 25 ) has certain drawbacks and 
enforces several assumptions for the analysis. This has mainly to do with the problem to 
determine P' and Z independently. The physical reason is that both lead to a reduction of 
the Andreev current and diminish the conductance change in G = dl/dV. The model fails 
to distinguish whether it is high P or high Z that causes the depression of conductance at 
small bias. This problem is evaded by applying a different theoretical approach. 26 
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Cuevas and coworkers 22,26 developed a model based on quasiclassical Green functions. 
The central quantities of the model are two transmission coefficients T-f-4 = |tt,J 2 - There- 
fore, we will refer to this model as the r^-r^-model throughout this paper. The transmission 
coefficients contain all microscopic properties relevant for the transport through the con- 
striction, i.e., they account for the majority- and minority-spin bands in the ferromagnet, 
the electronic structure of the superconductor, and the interface, and = a/1 — t^j,, 
respectively, are the spin-dependent transmission and reflection amplitudes, respectively, 
entering the normal-state scattering matrix S which supplements the boundary conditions 
of the theory. Of course, the restriction to a single conduction channel per spin direction is 
a rough simplification of the point-contact, but it is finally justified by the agreement with 
the experiment. 22 ' 23 Following the calculation by Cuevas and coworkers the spin- dependent 
current through the S/F point contact can be separated in two spin contributions, 

IsF = If + h 
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and each contribution can be written in its BTK form 

where f(E) is the Fermi function, A a (e) and B a (e) are the spin-dependent Andreev re- 
flection and normal reflection probabilities, respectively, and a =t or \.. A a (e) (B a (e)) is 
calculated from the spin-dependent transmission (reflection) amplitudes, and finally, the 
zero-temperature conductance of the S/F contact adopts the form 

, 2 ( in pV < A 

4e I (l+r t r 4 ) 2 -4r t r 4 ( e y/A) 2 C - a 



eV > A 



~~ h \ ¥i+(^r¥t)\/''(A/^) 2 

[(l-r t rJ+(l+r t r ; )Vl-(A/eV)^ 

while the normal-state conductance is given by 

e 2 

It is obvious that the Andreev spectra are determined by a set of three free parameters r^, 
rj,, and A. The current spin polarisation P in this model is defined by 

p = kt - n\ 
r t + n' 

and can be determined from the fit parameters of an experimental Andreev spectrum. We 
note, that this expression is symmetric with respect to and tj,, therefore, one cannot 



assign a transmission coefficient to the majority or minority charge carriers in the ferromag- 
net. However, we expect the high transmissive coefficient to correspond to the minority 
electrons. In the absence of spin polarisation, i.e., P = for a N/S contact with t-j- = tj,, 
above formulae reduce to the well-known BTK result 9 . 
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FIG. 1. Normalized differential conductance G/Gn vs eV/A. The curves are calculated with the 
T|-T|-model (see text) for different polarisations at T = 0.02T C . For clarity, the curves are shifted 
successively upwards by 0.2 units. 

Figure 1 displays a set of normalized conductance curves for T = 0.02T C , = 0.99, and 
r t = 0.99, 0.53, 0.42, 0.33, 0.25, and 0.17 which equals P = 0, 0.3, 0.4, 0.5, 0.6, 0.7 from top 
to bottom. The shape of each spectrum is unambiguously determined by a set of r-f-, rj.. For 
P = (top curve) the curve reproduces the well-known BTK result where Andreev reflection 
causes a doubling of the normal-state conductance for energies eV < A. We note that the 
characteristic double-peak feature at \eV\ = A caused by the reduction of the conductance at 
low bias originates from a small fraction of charge carriers ordinarily reflected at an interface 
barrier. In the case of two spin-dependent transport channels it is intuitive to consider the 
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high-transmittive spin channel to get a measure for the fraction of charge carriers ordinarily 
reflected at the interface barrier. From the transmission coefficients rj. = 1/(1 + Zf) one 
derives Z± = a/(1 — T^jJ^i- For the upper curve which corresponds to the unpolarised BTK 
case one gets Z^ — 0.1 as the parameter equivalent to the BTK interface parameter Z. 




eVM 



FIG. 2. Normalized differential conductance G/Gn vs eV/A. The curves are calculated with the 
T-f-Ti-model (see text) for the same polarisation P = 0.4 at T = 0.22T C . For clarity, the curves are 
shifted successively upwards by 0.4 units. 

While the shape of the spectra is unambiguously determined by a set of t-j-, t±, the spectra 
for same polarisation can look quite different. Fig. 2 shows four curves calculated for different 
sets of r-f, rj, which all result in P = 0.4. The calculations have been performed for a finite 
temperature T = 0.22T C . The high-transmission spin-channel seems to be decisive whether 
the curve shape appears more point- contact- like or more tunnelling-like. 

Before we compare curves calculated by both models, let us first check the validity of 
the fitting procedure. For this purpose we measured point-contact spectra of S/N point 
contacts in a 4 He cryostat down to T = 1.4 K and fitted them with both models. The 
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FIG. 3. Point-contact spectra of S/N contacts together with fits according to Martin- Rodero et 
al. 2e (dashed line) and Mazin et al. 25 (solid line). Upper panel: Pb/Cu contact at 1.4 K with 
R = 2.7 SI. Lower panel: Pb/Pt contact at 2.1 K with R = 7.5 CI. For the fit parameters see text. 



S 



PCs have been established in the edge-to-edge configuration with a sharpened Pb electrode 
and a normal-metal electrode made from Cu or Pt. Therefore, for both fits we expect 
P' = P — 0. Figs. 3a and b show the experimental data together with the fits according to 
the r^-T|-model 26 (dashed line) and the dispartment model 25 (solid line). In both cases we got 
almost perfect agreement with the data. For the Pb/Cu contact with R — 2.7 T2 measured 
at T = 1.4 K (upper panel) we obtained A = A(T = 0) = 1.38 meV, r t = t± = 0.814 
as parameters for the T|-T|-model, and A = 1.38 me V, Z = 0.475, P' = 0.011 for the 
dispartment model. For the Pb/Pt contact with R = 7.5 Q measured at T = 2.1 K (lower 
panel) the parameters are A = 1.48 meV, r/A = 0.10, r-j- = 0.811, t± = 0.819 for the r-j-- re- 
model, and A = 1.42 meV, r/A = 0.10, Z = 0.488, P' = 0.015 for the dispartment model. 
In both cases we found within the experimental error a good agreement of P ~ P' ~ 
as expected for these non-magnetic metals. The energy gap of Pb determined from the 
fits coincides fairly good to the gap value reported in literature 29 . We note that a small 
broadening of the Pb/Pt spectra caused by inelastic scattering in the contact region is 
accounted for by introducing the Dynes 30 parameter T which is of the order of 5-10% of Ao- 
In the next step we compare curves calculated with the simple dispartment model to 
those calculated with the Tf-T^-model for nominal same polarisation P = P'. Fig. 4 displays 
a set of normalized conductance curves calculated for T = 0.1 K, P = P' = 0, 0.2, 0.4, 
0.6, 0.8, and 1, and Z = 0.3 for the dispartment model, and rj. = 0.917 for the r-f- remodel, 
respectively, which corresponds to Zj. = a/(1 — 74.) /rj. = 0.3. In the extreme case P = P' = 
which describes a S/N contact the calculations perfectly agree with each other, as well as 
for the other extreme case P = P' = 1 which describes a halfmetallic F/S contact. For the 
latter, there is only a small difference in the vicinity of the coherence peaks at |eV/A| = 1. 
However, at intermediate values with increasing polarisation the Andreev signal is much 
faster suppressed for the dispartment model than for the r-f-r^-model. Our comparison 
discloses a notable difference of both quantities which makes questionable a contrasting 
juxtaposition of P and P 1 values derived from the analysis of experimental data by one of 
these models. 

In order to illustrate this difference on experimental data we used both models to fit 
the same set of data measured on Al/Fe nanostructured point contacts. 23 The PCs were 
fabricated by structuring a hole of 5 - 10 nm diameter with electron-beam lithography and 
subsequent reactive ion etching into a 50-nm thick Si3N 4 membrane, evaporating on one 
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FIG. 4. Comparison of Andreev spectra with P = 0, 0.2, 0.4, 0.6, 0.8, and 1 calculated according 
to the Tf-T^-model (black lines) and the dispartment model (red lines). For clarity, the curves are 
shifted successively downwards by 0.4 units. 

side a 200 nm Al layer, and on the other side a 12-nm thick Fe layer and a Cu layer of 
thickness <i Cu = 188 nm as a low-ohmic electrode. 31 Fig. 5 displays the normalized conduc- 
tance spectra of six different nanostructured PCs with contact resistances between 2.70 
and 24.2 measured in a dilution refrigerator at T = O.IK together with fitting curves 
calculated with the t-j-t^- model (left panel) and the dispartment model (right panel). Both 
models perfectly describe the data, minor deviations are observed only at |eV/A| ~ 1 where 
the experimental curves are more rounded probably caused by a leveling-off of the electron 
temperature due to heating by electromagnetic stray fields or a small pair-breaking effect by 
Fe. The corresponding fit parameters are listed in Table I. Within an uncertainty of 1% the 
same gap value Ao is found for both models, however, there is a notable difference in the Z 
parameter, which is a factor 2-3 higher in the dispartment model, and the spin polarisation 
P' which is lower. Although the origin of Z is not clear at all, in both models it subsumes 
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FIG. 5. Point-contact spectra of Al/Fe contacts at O.IK together with fits according to the r+- 
r^-model (left panel) and according to the dispartment model (right panel). 

all ordinary reflection of charge carriers that occurs at the interface for P = P' = 0, e.g., 
reflection caused by an insulating interface layer, lattice imperfections, Fermi velocity mis- 
match, etc.. For P and P' > the situation is less apparent. The physical reason is that Z 
and P' both lead to a reduction of the Andreev current. The dispartment model obviously 
fails to distinguish whether a high spin polarisation or a high barrier causes the depression 
whereas for the 7f-rj,-model as per definition only the conductance channel not affected by 
the suppression of Andreev reflection due to polarisation is considered to determine ordinary 
reflection. 

Another important aspect is that the previously reported dependence of the spin polar- 
isation on the contact size 23 is robust against the model used to extract P. Independently 
of the model there is a clear reduction of the spin polarisation with decreasing contact re- 
sistance Rn, i- e. increasing contact radius a. The reduction of P has been allocated as 
being due to spin-orbit scattering in the contact region with a constant scattering length £ so 
modelled by a simple exponential decay 23 , P(a) = Po ex P ( — a /4o)- A spin-orbit scattering 
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Sample 


Rn 


A 


7"t T i Z i P 


A 


Z P' 


No. 


(ft) 


(meV) 




(meV) 




1 


2.68 


0.174 


0.371 0.983 0.132 0.452 


0.176 


0.337 0.396 


2 


6.98 


0.175 


0.362 0.984 0.128 0.460 


0.176 


0.338 0.407 


3 


7.29 


0.157 


0.349 0.993 0.083 0.480 


0.158 


0.272 0.444 


4 


9.59 


0.190 


0.361 0.984 0.128 0.463 


0.191 


0.342 0.407 


5 


18.4 


0.166 


0.348 0.997 0.054 0.482 


0.168 


0.218 0.497 


6 


24.2 


0.174 


0.343 0.994 0.078 0.487 


0.174 


0.262 0.494 



TABLE I. Fit parameters Ao, t+, and t\ of the r^-r^-model and Ag, Z, and P' of the dispartment 
model for conductance spectra of nanostructured Al/Fe point contacts. From the transmission 
coefficients r-t- and Ti the interface barrier Z^ and the current spin polarisation P have been calcu- 
lated. 

length i so = 255 nm has been obtained for the analysis with the r^-r^-model. The same 
systematic trend of P'{a) is found for the dispartment model albeit with a lower value for 
£ so . For small contacts both models result in almost the same spin polarisation. 

In summary, we have discussed possible reasons for the scatter of polarisation values found 
for the spin polarisation measured by Andreev reflection in point-contact experiments. We 
showed that the scatter is partially caused by the models used to extract the spin polarisation 
from the data, and partially caused by intrinsic mechanisms like the spin-orbit scattering in 
the contact region. 
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